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Abstract 



We analyze the spontaneous baryogenesis and charge transport mechanisms sug- 
gested by Cohen, Kaplan and Nelson for baryon asymmetry generation in extended 
versions of electroweak theory. We find that accounting for non-perturbative chirality- 
breaking transitions due to strong sphalerons reduces the baryonic asymmetry by the 
factor (mt/nT) 2 or aw, provided those processes are in thermal equilibrium. 
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The popular mechanisms for electroweak baryogenesis discussed in the literature 
are the spontaneous baryogenesis mechanism and the charge transport mechanism, 
both suggested by Cohen, Kaplan and Nelson |], |3|. It has been advocated in 
fl], 0, § that these mechanisms produce an asymmetry parametrically larger that the 
mechanism associated with the decay of topologically non-trivial fluctuations of the 
gauge and Higgs fields during the electroweak phase transition, considered in refs. ||, 
3 (for earlier suggestions, see refs. @, 0, §])• In the latter case, the asymmetry produced 
is proportional to the square of the ratio of the top-quark mass to the temperature 
while for the first two mechanisms it was argued that this suppression is absent 
|1|, ||, [|. The key observation made in [p], ||, |3| is that in the presence of a non- 
zero fermionic hypercharge density, py, the equilibrium value of the baryonic charge, 
generally speaking, is not equal to zero: 

(B) = Apy, (1) 

where A is a number determined by the particle content of the theory and A ^ in 
the limit m t — > 0. 

In this paper, we show that under some circumstances A may be zero in the massless 
quark limit. Moreover, we will argue that this is most probably the case for the two 
Higgs doublet theory, considered in |TJ, 0. 

In this model, in addition to the KM CP-violation, there is an explicit CP-violation 
in the Higgs sector. The vevs of the Higgs fields at zero temperature have the form: 

(4 >1 ) = (0,v 1 e ie ), (&) = ((), -ya), (2) 

where 9 is the relative phase, which cannot be rotated away by a gauge transformation. 
Due to explicit CP-violation in the Higgs potential this phase is space-dependent inside 
bubble walls forming at the electroweak phase transition (it would be time-dependent 
for a spinodial decomposition phase transition). Since quarks are getting their masses 
from interactions with the Higgs fields, the most important interaction is the top-quark 
coupling^, 

Ly = ftQsUsfa. (3) 

Here Qi are the left fermionic doublets, C/j, Di are right quark fields, and i is the 
generation index. The bottom-quark couplings can be safely ignored provided the 
thickness of the domain wall Iq or the duration of the spinodial decomposition phase 
transition tq is smaller than the rate of chirality-changing transitions associated with 
the bottom Yukawa coupling; that is, if 

ID 4 

l < (aj'T)- 1 ~ — , (4) 

where T is the temperature. This is always the case for current estimates of the domain 
wall thickness, ranging from 40/T H] to 1/T [TO. The same statement is obviously 
true for lighter quarks. 

3 We assume that the first Higgs field couples only to up-quarks, while the second is coupled only 
to down-quarks. This assumption is, however, inessential in what follows. 
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We will clarify our point first within the context of the spontaneous baryogenesis 
mechanism. Let us take for simplicity the spinodial decomposition phase transition. 
As in PI, we assume that the rate of chirality- flipping transitions for the top quark is 
larger than the typical inverse-time scale Tq" 1 of the variation of the phase 9 during this 
transition. Now, following we make a hypercharge rotation of the fermionic fields 
in such a way that time-dependence disappears from the Yukawa coupling (|]). This 
change of variables introduces the following modification of the effective lagrangian: 



L eff = L + 6Y F 
where Yp is a fermionic hypercharge operator, 

3 



8=1 



1 - 4 - 2 - 

-QaoQi + qUiioUi - -A70A 



La Li - 2Eij Ei 



(5) 



(6) 



We want to check whether in thermal equilibrium with respect to fermionic number 
non-conserving processes in a theory with effective action (||) the baryonic number is 
nonzero. To determine the number A defined in ([!]), the following standard procedure 
must be used (see e.g. ref. |ITfl ). One has to define the complete set of conserved 
charges X« and construct the most general equilibrium density matrix with the help 
of those charges, introducing chemical potentials //j for each of them. The partition 
function is 

1 (H-6Y F -J2^X; 



Z = Tr exp 
while the density matrix is 

1 

p = — exp 



T 



(7) 



Now, since X, are conserved operators, their average must be equal to zero. This 
requirement fixes the chemical potentials /ij through the system of equations 



d_ 
Ik 



0. 



Then the baryonic number is just 



(B) = Tr[Bp}. 



(9) 



(10) 



Let us define the set of conserved numbers. It is sufficient to consider only SU(3) 
gauge singlets. We start from purely fermionic currents. In a model with massless neu- 
trinos (no right-handed neutrinos!), the total number of different fermionic currents 
is 15, the number of fermionic degrees of freedom. They are: 3 left quark currents, 
QiJtiQi, 6 right quark currents, IJcipUi and Dij^D^, 3 left leptonic currents, Lij^Li 
(Li is the left leptonic doublet), and 3 right leptonic currents, E^^Ei. Not all of these 
currents are conserved. One has to take into account the following processes: 
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(i) Perturbative chirality-changing transitions due to Yukawa interactions (Q). The 
rate of these processes is estimated to be r Y ~ 30/T 0. This decreases the number 
of conserved currents by one: instead of the pair Q^^Qz and U^^U^, we have a con- 
served linear combination Q^^Qz + U^j^U^. 

(ii) Non-perturbative chirality-breaking transitions due to strong interactions. It 
is well known that the quark axial vector current has an anomaly and therefore is 
not conserved. The rate of chirality non-conservation at high temperatures T str0 ng is 
connected with the rate of topological transitions in QCD ("strong" sphalerons, JI2]]), 

9Q 5 _ 12-6 

~ rps 1 strong^ 5; K 11 ) 

where Q$ is the axial charge. The factor of 12 comes from the total number of quark 
chirality states, the factor of 6 from the relation between the asymmetry in quark 
number density and the chemical potential, 

:\n F (e,p)-n F (e,-ii)] = -—(l-—j=), (12) 



(2tt) 3L v v ' ^ n 6 v 2tt 2 T 2 

where n F {e, /i) = l/(exp(^r £ ) + 1) is the Fermi distribution, e 2 = k 2 + m 2 . The rate 
of the strong sphaleron transitions is related to the rate of weak sphaleron transitions 
through 

^strong = ~( — ~ sph = ^ K ( a sT) 4 , (13) 
6 Otw 6 

where k is the usual parameter characterizing the strength of the electroweak sphaleron 
transitions in the unbroken phase. The characteristic time of these transitions is there- 
fore | 

Tstron9 = «' (M) 



Using the conservative lower bound on k derived in lattice simulations |13|], k > 0.5 
(see also discussion in ref. |H|), we obtain a conservative bound T strong < 100/T. With 



k ~ 20, so estimated by a different method in ref. fToH , we obtain T strong ~ 2.5/T. 
Therefore, the rate of strong sphaleron transitions is comparable to or even larger than 
the rate of chirality-flip transitions through the Yukawa coupling of the top quark. 
Hence, these processes must be taken into account. This decreases the number of con- 
served currents by oneFl 



(iii) Anomalous fermion number non-conservation must be taken into account |T7 
It decreases the number of conserved currents by one. 



4 Note that influence of strong sphalerons on a different mechanism for electroweak baryogenesis 
has been discussed also in ref. [|l6| . 
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Therefore, a complete set of conserved anomaly-free fermionic currents can be rep- 
resented as 



Qil^Qi + R - L L , i = 1, 2, 


(15) 


Qzl^Qz + U 3 ^U 3 + R/2 - L L , 


(16) 


Lij^Li - L L /3,i = 1,2, 


(17) 


Ea^E^i = 1,2,3, 


(18) 


A7„A-fl/2,i = l,2,3, 


(19) 


Un^Ui - u 2lfl u 2 , 


(20) 



where 

3 2 

L L = ^LijpLi, R = Y,UilJJi. (21) 

i=i i=i 

In addition to the quantum numbers associated with these currents, there exists a 
conserved operator containing scalar fields. This is the familiar hypercharge 

Y = Y F + Y S , (22) 

where 

Y s = -i[<j>\V Q <k - (Z>o0i) Vi] - i[l - 2] (23) 

is the scalar field contribution. Finally, one should add the third component of the weak 
isospin operator T3 . As in [0], [| , we assume that the classical motion of the scalar field 
does not introduce any non-zero hypercharge or T 3 density. Therefore the hypercharge 
density is equal to zero for a uniform spinodial decomposition phase transition. 
Now one can check that, in the massless quark approximation, eqs. (0, ^|, |9], 



UD imply that the equilibrium value of the baryonic charge is zero, (B) = 0. In order 
to get this result, we used, following ref. pjjjl , the particle spectrum of the unbroken 



phase rather than the broken one. 

The reader may wonder why we did not use the electric charge operator and the 
physical particle spectrum of the broken phase instead of Y and T 3 []. The reason 
is that the use of the physical spectrum in such a calculation would correspond to 
a computation of the free energy of the system in the broken phase in the one-loop 



approximation in a unitary gauge. It is known (T^, ^(| |2]J, that perturbation theory in 
a unitary gauge is a very delicate thing (e.g., it does not converge 01). In particular, 
it gives an incorrect result for the critical temperature of the phase transitions in gauge 
theories. 

It is clear that for sufficiently small vacuum expectation values of the Higgs field, all 
effects associated with spontaneous symmetry breaking are just the mass corrections 
which can be neglected in the high-temperature approximation in any renormalizable 
gauge. In the worst case, the expansion parameter may be (v/T) 2 . In our case, the 
vacuum expectation value of the Higgs field v is bounded from above by the requirement 



5 Actually, the zero result is reproduced also with Q instead of T3 and Y and the particle spectrum 
of the broken phase. 
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that the sphaleron processes are sufficiently fast, so that v < gwT. Therefore, the use 
of the particle spectrum of the unbroken phase is perfectly justified. 

The zero result is also reproduced when fermionic number non-conservation is out 
of thermal equilibrium. Then one can write an equation for the baryonic number 
evolution [|T8| (see also refs. |1|, [2|. [3|] ) 



dt T dB T 1 ' 

where F(B) is the free energy of the system characterized by zero values of all conserved 
charges but non-zero baryonic density B, T sp h is the rate of sphaleron transitions, and 
N = 3 is the number of fermionic generations. In our case, 

dF 

q^ = ^b-B, (25) 

so that B stays zero all the time. 

We stress that the inclusion of strong sphalerons was essential for this result. Strong 
sphalerons have the physical effect of mantaining the same chemical potential for left- 
and right-handed baryonic numbers. If instead one neglects these processes, then the 
system of charges (^) has to be supplied with an extra charge, say R, while the density 
matrix will contain an additional chemical potential. One can check that if one chooses 
this set of conserved charges, the result is indeed non-zero. Which solution has to be 
used? The answer depends on the ratio of the typical time scales. If r <C r strong , then 
strong sphalerons are irrelevant, and the asymmetry indeed does not contain Yukawa 
coupling suppression provided that t y is smaller than the typical time scale r of the 
Higgs phase change in the spontaneous baryogenesis mechanism. If, on the contrary, 
to > t strong or t ~ t 'strong, then the asymmetry vanishes in the massless approximation. 

In the original calculation of Cohen, Kaplan and Nelson fl], [| a non-zero result for 
the baryonic asymmetry was obtained, since strong sphaleron transitions have been 
neglected and only B — L, Q, and B 3 — §(-£>i + B 2 ) have been included as conserved 
charges. However, in presence of strong sphalerons transitions, their solution corre- 
sponds to non-zero averages for some of the approximately conserved charges defined 
ineqs. (|L|-|20|). 

Is the conclusion that (B) = fatal for the spontaneous baryogenesis mechanism? 
In fact, this is not the case when quark mass corrections are included. Mass corrections 
can be taken into account in perfect analogy with the case of the leptonic asymmetries 
discussed in refs. |22|, |18| (for a later discussion see ref. p3|). To make the computation 



simpler (and to make clearer why we previously obtained zero), we observe that in fact 
the number of independent chemical potentials (14) can be reduced to a smaller number 
(6) with the use of flavour symmetry. Namely, our lagrangian has quite a large global 
symmetry group SU{2)q x SU(3)d x SU(2)u x SU(3)l x SU(3)e, where the group 
labels denote the fields on which the corresponding group is acting. The conserved 
charges invariant under this symmetry are 

2 

Ai = Y, QHnQi + 2R ~ 2L l, (26) 
1=1 
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A 2 = Qzl^Qs + U 3 yJJ 3 + R/2 - L L , (27) 

3 

^3 = E Day. A - 3R/2, (28) 

3 

A^^Ea^, (29) 
i=i 

hypercharge Y and weak isospin T 3 . We denote the corresponding chemical potentials 
by Hi, i = 1, ... , 4, /iy and /iy. They are to be found from the equations 

(A) = (Y) = (T 3 ) = 0. (30) 

Let us first solve this system neglecting the mass of the top quark. In order to show 
that the baryonic density is zero, one can consider only the total baryonic density and 
the left-handed leptonic density. The total baryonic density is 

(B) = 4[{9 + fi Y ) + 2^ + (31) 

o 

while the total left leptonic density is 

T 2 

(L L ) = -60 + ^+2//! + /^—. (32) 

o 

Then from (B — Ll) = 0, it follows (B) = 0. For future reference, we give here the 
complete solution of eq. (30) in the massless approximation: 

Mi = ~2i^ + Py), ^2 = + Py), Ms = ^0 + Hy), (33) 

14n 

^ = 2(6 + fi Y ), fi T = 0, (6 + fiy) = 9 + 14 B W8 , (34) 

where n s is the number of scalar doublets. 

The exact proportionality between baryonic number and left-handed leptonic number 
is, however, an artifact of the massless approximation. Using eq. fll2j), the baryonic 
density becomes 

/R\ 9Us 2a /orN 

\B) = tt^Py = ; 6, 35 

v 1 20vr 2 T 2 ^ 10(9 + Un s )7i 2 V ; 

where py is the fermionic hypercharge density. This is our final result for the equilib- 
rium density of the baryonic number in the background of a scalar field with slowly 
changing phase. 

If strong sphalerons were not in equilibrium, we should add a new approximately 
conserved charge, for instance A 5 = R, to our previously defined set of charges, see 
eqs.(26)-(29). Following the same procedure used to obtain eq.(35), we would find: 

(B) = TH. (36) 

6 + Wn s 
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Therefore the presence of strong sphalerons leads to a suppression factor y§§(^|0 2 , 
where we have taken n s = 2; this corresponds numerically to a suppression of about 
3 • 10" 2 for f t ~ 1, v{T) ~ g w T. 

The approximation that anomalous baryon number violation is in equilibrium is 
unlikely to be correct. Thus, following ref. [|1], we consider the evolution equation for 
the baryon asymmetry density, eq. (24). We can compute /!# using the same procedure 
followed above, adding B to the set of conserved charges and imposing, besides eq. (30), 
the further constraint (B) = 0. We find, if strong sphalerons are in equilibrium, 

^ " 4(9 + 14n s ) vr^' (3?) 
and, if strong sphalerons are not in equilibrium, 

= -irfkf (38) 

The mass suppression present when strong sphalerons are in equilibrium amounts to a 
factor f§(^) 2 ~ 2- 10" 2 . 

We can now consider the case of baryogenesis with the charge transport mechanism. 
As shown in ref. 0, the CP- violating interactions of the fermions with the moving thin 
bubble wall lead to a non- vanishing flux of particles carrying some quantum number X. 
If, as in the original work of ref. [[J , we identify X with the hypercharge, then it is easy 
to realize that (B) = 0, when strong sphalerons are in equilibrium and all fermions are 
massless. The calculation is analogous to the one for spontaneous baryogenesis, with 
the only difference that, instead of the source 9 we have to consider a non-vanishing 



Y density. However, as shown by Khlebnikov [p4 |, the Debye screening prevents the 



transport of the gauged charge Y over distances larger than about 2/T. 



Therefore the recipe proposed in ref. [25] is to identify as the transported charges X 



only global charges which are orthogonal to gauge charges. We then define 

B' = B + a B Y, A^ = Ai + OiY, i = l,...,4, (39) 

with as and ccj chosen such that B' and A[ are orthogonal to Y, i.e., in a massless 
quark approximation: 

112 2 1 

-a B = — ai = ~a 2 = -a 3 = -a 4 = . (40) 

By introducing chemical potentials for the charges (26)-(29), Y, and B, we can write 
the asymmetry density for each species of particles as: 

P=Y (otPy + Qb^b + lAil^ij ■ (41) 
The chemical potentials /iy, //g, and //j are then computed from the equations 

(Y)=0, ^ = ^ = r^0, (42) 
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which correspond to a situation in which the gauged charge Y is screened, but the Y 
components of the global charges can carry a non- vanishing density n x . The solution of 
eq. ( pE2"D yields //# = and therefore anomalous baryon number violating interactions 
are unable to generate any baryon asymmetry. One can also check that 7^ 0, if 
strong sphalerons are not in equilibrium. 

The cancellation in the case of charge transport baryogenesis may seem more severe 
than in the case of spontaneous baryogenesis since the relevant processes occur in 
the unbroken phase, where the Higgs field condensate is absent. There are two types 
of corrections which save the situation. Just as in the leptonic case considered in 
Yukawa and gauge radiative corrections are important. The equation for the 



asymmetry number density of each chiral fermion in the unbroken phase, including 
radiative corrections, can still be written as in eq.(12), with the replacement: 



m 2 



C s 9s + ( c w + —r sin 2 9 w )9w + tt 



T 



2 



(43) 



where C s is I for quarks and for leptons, Cyy is | for weak doublets and for 
singlets, Y is the hypercharge quantum number, and / is the Yukawa coupling, taken 
to be non-vanishing only for the top quark. The solution of eq. fl42|) now gives: 

which, as noticed above, vanishes in the limit of negligible radiative corrections. The 
result of eq. (^) would correspond to a reduction of the baryonic asymmetry in 
comparison with [FJ, |3| by the factor ~ 12871-2 (9+i4n ) ft ~ 3 • 10 ^, if other effects are 
absent. These more important effects are associated with electroweak corrections which 
break the degeneracy between left and right top quarks in the unbroken phase resulting 



in a net baryonic current Jb (see refs. p6], [141 ) related to the hypercharge current Jy 
through 

T % PL -PR T ,._v 

Jb = ■ Jy, (45) 

3 Pl + Pr 

where pi and pn are the momenta of incident left and right quarks corresponding to 
the same energy. Now, pi and Pr are different since left quarks interact with SU(2) 
gauge fields while right quarks do not. With the use of eq. fl4~3|) we get a larger baryonic 
chemical potential, 

where p > m t is a typical transverse momentum of the quarks contributing to charge 
transport. So the final estimate of the suppression of the asymmetry in this case 
is a factor of about 10~ 2 . Therefore, the charge transport mechanism (as well as 
the "topological" mechanism can still explain the baryonic asymmetry since, 

according to the estimate of ref. 0, ub/s can be as large as 10~ 7 . With strong 
sphalerons taken into account, this number should be converted to 10~ 9 , which is still 
larger than observation. 



^ —nr— ' ( 46 ) 



s 



To conclude, we have shown that the baryonic asymmetry in the two Higgs doublet 
model produced by any of the mechanisms considered in |l], H, H contains a parametric 
suppression associated with the top Yukawa coupling or electroweak coupling due to 
the existence of strong sphalerons. Strong sphalerons have not been taken into account 
in a number of papers on electroweak baryogenesis and may change their conclusions. 

One of us (M.S.) is indebted to the High Energy Physics Group at Rutgers University, 
where part of this this work has been done, for kind hospitality. We thank G. Farrar, 
D. Kaplan, S. Khlebnikov, A. Nelson and N. Turok for useful discussions. 
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